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Abstract: A three-dimensional unsteady reactive Navier-Stokes code is developed to 
study the ignition and flame spread above liquid fuels initially below the flashpoint tem- 
perature. Opposed air flow to the flame spread due to forced and/or natural convection 
is considered. Pools of finite width and length are studied in air channels of prescribed 
height and width. Three-dimensional effects of the flame front near the edge of the pool 
are captured in the computation. The formation of a recirculation zone in the gas phase 
similar to that found in two-dimensional calculations is also present in the three-dimensional 
calculations. Both uniform spread and pulsating spread modes are found in the calculated 
results. 

Introduction: Schiller, Ross, and Sirignano [1] studied ignition and flame spread above 
liquid fuels initially below the flashpoint temperature by using a two-dimensional computa- 
tional fluid dynamics code that solves the coupled equations of both the gas and the liquid 
phases. Their computational studies and analysis identified the mechanisms for uniform and 
pulsating flame spread and studied the effects of gravity level, pool depth, fluid properties, 
and chemical kinetic coefficients on flame spread across liquid fuel pools. Pulsating flame 
spread was attributed to the establishment of a gas-phase recirculation cell that forms just 
ahead of the flame leading edge because of the opposing effect of buoyancy-driven flow in the 
gas phase and the thermocapillary-driven flow in the liquid phase. Schiller and Sirignano [2] 
extended the same study to include flame spread with forced opposed flow in the gas phase. 
A transition flow velocity was found above which an originally uniform spreading flame be- 
comes pulsating. The same type of gas-phase recirulation cell caused by the combination of 
forced opposed flow, buoyancy-driven flow, and thermocapillary-driven concurrent flow was 
found to be responsible for the pulsating flame spread. Ross and Miller [4] and Miller and 
Ross [5] performed experimental work that generally corroborate many of the computational 
findings in References 1 and 2. 

In this paper, we extend previous two-dimensional and axisymmetric studies in References 
1, 2, and 3 to three-dimensions. Figure 1 shows the geometry of the three-dimensional model. 
The mathematical formulation of the problem and assumptions are the same as those in 
References 2. 
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Results and Discussion: Four cases of different initial temperature T 0 at normal 
gravity and zero gravity are studied: (1) T 0 = 21°C at normal gravity; (2) T 0 = 21°C at zero 
gravity; (3) T 0 = 25°C at zero gravity. (3) T 0 = 15°C at normal gravity; The flow velocity 
at the exit of the wind tunnel is fixed at 30 cm/s for all four cases. 

Figure 2 shows the contours of fuel consumption rate 0.89 seconds after ignition in the 
x—y plane at zjL — 0.02 for Case 1. Figure 3 shows contours of temperature and flow velocity 
vectors in the same plane. The spatial coordinates presented here are non-dimensionalized 
by the total length of the wind-tunnel, i.e., 36 cm. Near the surface of the liquid pool, surface 
tension pulls the fluid forward ahead of the flame while the forced flow’ velocity goes against 
this motion. This causes small vortex cell in front of the flame as shown by the velocity 
vectors shown in Figure 3. A top view 7 of the flame ( x — z ) plane at a dimensionless distance 
of y/L = 0.0022 from the liquid fuel surface is shown in Figure 4. The flame at the corner 
curves around and bends at the edge of the fuel pool. The flame reaches out of the fuel pool 
and attaches to the wind-tunnel floor by the pool edge. Figure 5 shows the contours of fuel 
consumption rate in a cross-section at x/L = 0.805. Figure 6 plots the velocity vectors and 
contours of temperature in the same plane. Buoyancy effects pull the flow into the flame 
from the side very near the wind tunnel floor and close to the flame whereas the flow 7 expands 
out from the pool away from t he wind-tunnel wall. - 

Figure 7 show’s the calculated history of the position of the maximum fuel consumption 
rate location for all four cases within the computed time period. At the time of the waiting 
of this abstract, computations for Cases (2) and (4) have been performed for the longest 
time to show 7 pulsating features of the flame spread. 

In summary, A three-dimensional computational code has been developed to solve the 
coupled reactive flow 7 problem of forced gas phase flow 7 over liquid fuel. Three-dimensional 
effects of the flame front near the edge of the pool are captured in the computation. The flame 
curves backwards near the pool edge and extends beyond the edge of fuel pool to be attached 
to the wind-tunnel wall. Thermal expansion and buoyancy effects c reate a s tr eamwise vortex 
near the two edges of the liquid pool. Despite the three-dimensional effects, computations 
show 7 that the basic features oT the flame propagation process is the same as tliose in the two- 
dimensional simulations. The formation of a recirculation zone in the gas phase is necessary 
for pulsating spread of the flame. Further studies w’ill include more detailed simulation of 
the propagation of the flame over the w'hole pool surface and also variation of parameters 
such as forced flow speed, initial temperature, liquid pool width and depth. 
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Figure 1: Geometry for the three-dimensional liquid fuel pool in a wind-tunnel with forced 
flow 
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